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Abstract 

We compute holographically the diffusion constant of supercharges in N=4 SYM 
at finite chemical potential for the R-charge, by solving the equations of motion for 
the transverse mode of the gravitino in the STU black hole in 5 dimensions. We 
consider the case of one charge and three charges, and we present analytical solutions 
for small values of the charges and numerical solutions for arbitrary values. We 
compare our results with other known results in 4 dimensions. 
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1 Introduction 

A large body of literature has been devoted to the study of the transport and 
hydrodynamic properties of strongly coupled gauge theories using the AdS/CFT 
correspondence, starting from [HE] [3]; see [I] for a recent review. 

One of the most striking results of these investigations has been the discovery 
that the ratio of the shear viscosity to the entropy density is a universal quantity, 
that has the same value of 1/Att in every theory that admits a description in terms 
of Einstein gravity coupled to matter fields; the universality is spoiled if the La- 
grangian has higher-derivative terms, which means generically, given the dictionary 
of the holographic correspondence, that it holds only in the limit of infinitely strong 
coupling of the dual field theory. See [5] for a review. The underlying reason is the 
uniqueness of Einstein gravity, or in other words the universal coupling of gravity 
to all matter fields; as a consequence, it can be shown that in any planar black hole 
background, the transverse graviton (i.e. with polarization parallel to the direc- 
tions in which the horizon extends) satisfies the same equation of motion; the shear 
viscosity is related to the absorption of this mode by the black hole horizon. 

This result was put in a more general framework by |6j , who showed that for some 
transport coefficients, the relevant correlator has no radial flow, or it is independent 
of the radial direction in the low-energy limit; this happens whenever the fluctuating 
mode that couples to the operator is a massless scalar. In that case, the transport 
coefficient that is in principle computed on the boundary can also be computed on 
the horizon, and so typically it has less dependence on all the details of the model; 
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but it is not enough to imply universality as there can still be a residual dependence 
of the "effective coupling" at the horizon. 

Most of the work on the subject has been devoted to bosonic correlators, with rel- 
atively little attention paid to the fermions. In fact, it is natural to wonder whether 
the constraints imposed by supersymmetry may be sufficient to imply universality 
for the transport of supersymmetry charges. With this question in mind, one of the 
authors (GP) has computed the supercurrent diffusion constant in N=4 SYM |7J. 
In that paper the result was obtained by looking at poles of the retarded correlator 
of the supercharge density ; the imaginary part of the pole gives the attenuation 
coefficient which can be related to the diffusion constant by general hydrodynamic 
arguments. While this was the simplest way to arrive at the result, it was not opti- 
mal because it used the longitudinal channel instead of the transverse one; even for 
the graviton the universality would not be manifest in this channel. 

In the present paper we reanalyze the question by looking at the transverse mode 
of the gravitino. First we consider an uncharged gravitino in a generic black-brane 
background; the technique of radial flow allows us to derive a general formula for 
the diffusion constant; this formula reproduces the known result in the case of the 
AdS black hole; however even in this simple case the fluctuation is not massless in 
general, so we cannot read the value in terms of horizon data. 

Next we consider what happens when the gravitino is charged, as it is the case in 
gauged supergravity. The equations of motion in this case are much more complex 
and the radial flow is not useful. Therefore we consider only one particular case: 
the black hole in the so-called STU model. This is a solution of N=2 gauged 
supergravity, but it can be embedded in N=8 and as such corresponds to taking 
N = 4 SYM at finite chemical potential for the R-charge. The STU black hole 
can be charged under three commuting U(l) gauge fields. We consider the special 
cases of one charge and three equal charges. In each case we solve numerically the 
equations of motion and compute the diffusion constant as a function of the ratio 
of temperature over chemical potential. 

This is the plan of the paper: in section 2.1 we review the properties of the 
supersymmetry correlators in the hydrdynamic regime, and the setup for the holo- 
graphic computation; in section |2.3| we derive the diffusion constant for the un- 
charged gravitino; in section 3.1 we review the solution for the STU black hole and 
their thermodynamic properties; in sections |3.2 3.3| we detail the computation of 
the diffusion constant and give the results (in fig. 1-4). The appendices contain 
some details on the derivation of the Kubo formula for the supercurrents, and the 
normalization of the boundary action. 
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2 Supercharge diffusion at zero density 



2.1 Supercurrent correlator 

A supersymmetric field theory at finite temperature has hydrodynamic excitations 
corresponding to fluctuations of the supercharge density (see [8] and [9 J for more de- 
tails on hydrodynamics with supercharges) . In particular there is a supersymmetric 
sound wave that propagates with a speed and attenuation that are different from 
those of ordinary sound. The attenuation rate is governed by the diffusion constant 
of the supercharge; this has been computed in the N=4 SYM theory for the first 
time in [7] using holography. In this section we compute it in a different way, not 
using the sound attenuation but the transverse part of the correlator as explained 
below; along the way we present the definitions and techniques that will be used in 
the next section. 

The diffusion constant can be extracted by means of a Kubo formula (see ap- 
pendix |A|) from the retarded correlator of the supercurrents Sf, defined as: 



Gjf(fc) = J d±xe- ik -*ie{x G ){{Sf{x), £>)}> (2.1) 

Conservation of the supercurrents and superconformal invariance imply: 

k'Gij = fGy = drf = (2.2) 

The correlator of the supercurrents can be build out of projectors that implement 
these constraints. The projector on the transverse gamma-traceless part of a vector- 
spinor is 

1 ( $i \ ~ 1 



H = % - 3 (jK - 7 J " 3^(4/0, - Jd)k 3 (2-3) 

The correlator can be written as = P^M kl P\. At zero temperature, Lorentz 
invariance dictates the form of M to be M k i = A(k 2 ) f.r\ k i. At finite temperature 
or density Lorentz invariance is broken and the correlator can depend apart from 
k also on the velocity of the fluid u l = (1, 0, 0, 0). One can write another projector 
that is transverse both to k and u: 

Pli =P?i=0 

P jk = Sjk - 2 ( 7j " Ik ~ ~ 

where k l = {uj — /x,q). We define the longitudinal projector as P L = P — P T . The 
correlator can have three possible structures: P L P L , P T P T and P L P T + P T P L . In 
the following we choose q = (0, 0, q). With this choice the mixed correlator vanishes 
identically. The longitudinal and the transverse part of the correlator can be written 
as: 

G% = (P L )fM km (P L )mj 

G% = {P T )*M km (P T U ' 
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and the most general form of M and M allowed by the symmetries is: 

(2.6) 

There are additional relations between the coefficients that come from requiring ro- 
tational invariance at zero momentum. For the following we will need these relations 
between the components of the correlator in the low momentum limit: 

G22 = ^33 = — a '/2 £44 = 

G22 = G33 = -«'/6 G L U = -2o'/3 

The longitudinal part of the correlator was calculated in [7] . In the present article 
we focus on the transverse correlator. 



2.2 Transverse gravitino in blackbrane background 

For the holographic computation we have to consider a gravitino propagating in the 
near-extremal 3-brane background (see [TU] for the analogous computation at zero 
temperature). The near horizon geometry has the following metric: 

ir 2 T 2 R 2 R 2 
ds 2 = (-f(u)dt 2 + dx 2 + dy 2 + dz 2 ) + - - . du 2 (2.8) 

u 4:J{U)U Z 

where f(u) = 1 — u 2 . The boundary is located at u = and the horizon at u = 1. 
The bulk action for the gravitino is: 

S = j d^x^^^D^p - m^r^tf*,) (2.9) 



The covariant derivative acts on a spinor as: 

1 

4 

and the equation of motion is the Rarita-Schwinger equation: 



D^ = d^ + -ujf lab (2.10) 



T^Dy&p - mr""^ = (2.11) 
In the gauge ^ = the equations for the spatial components of the field read 

7 + ^TV^J { 7J 1 dt + l dj ) ^ k + ~^T 7 * k " ^ 7fc ^ 5 + W7^* = ° 

(2.12) 

where 4> a = , j a = e^T^ and i/i' = dip/du. 

Define the Fourier transform as: 

%(u,x) = ^ J dtke^y^k) (2.13) 



with k» = 2vrT(w,0,0,g). 

The field 77 = j 2/ tp2 — 7 3? / ; 3 has spin 3/2 under the 0(2) rotational symmetry that 
preserves k^ and therefore it decouples from the other components that have spin 
1/2. It is easy to see that this is the field that is dual to the transverse part of the 
correlator. Multiplying (2.12) for k = 2 by 72 and for k = 3 by 73 and subtracting 
we arrive at the equation satisfied by the transverse component of the gravitino: 

i 2 - 2 , 3 



>1 



-7 



+ 



uf 2uf 



-7 



f] 



(2.14) 



with ft = —iuj/y/f^ 1 + ic/7 4 . 

We cannot solve equation (2.14) analytically but we can solve it as an expan- 
sion in momenta. A useful remark is that equation (2.14) commutes with 723 = 
diag(l, —1, 1, —1) and therefore we can project on its eigenvectors. In the eigenspace 
723 = 1 (resp. —1) equation (2.14) becomes a matrix differential equation for (771, 7/3) 
(resp. (772,774))- 

We find the following solutions near the horizon : 

(1 - u)~Hf (-i, 1) 

(l-lt)-2 +i f(i,l) 

We have to impose incoming boundary conditions at the horizon in order to cal- 



(2.15) 



culate the retarded correlator pJJ.We can keep the incoming waves at the horizon 
by imposing the relation 77 + = —ir]~. For the two other components (772,774) the 
equations are the same after the exchange q —¥ —q. Repeating the analysis for these 
components we find the solution up to first order in momenta. The near boundary 
expansion of this solution reads: 



.•)u +«V*££ / 3( w , g ) 



/ ^ 4 H3 3 ^)( g -. 

u VH3 3 / 4 (u, + q) 1 + u 7 / 4 ^5(u>, q ) 
u 1 /^ 3 / 4 ^, q) + V2(q - (1 - V2l)u)a 
V u 1 /^ 3 / 4 ^, g) - V2(q + (1 - v/2£M7] / 



(2.16) 



where t = log(l + V2). We may identify the source of the field <p with the negative 
chirality part 77" = u 1 ^^ and identify \ an d X as the coefficients of the positive 
chirality expansion 77 + = u 3 ^x + u 7 ^X : 

( 3 3 / 4 [if3(uj, q)) + V2(q - (1 - v^)w)a] 
V 3 3 ^[i5(oo, q) - V2(q + (1 - v^Vh] 



(2.17) 



X 




33/4 6(u,q) 



X 



z3 3 / 4 (lj -q)a 
73 3 / 4 (w + q)-y 



(2.18) 



Multiplying <f> by $ = —i2irT diag(a; — q, uj + q) and neglecting terms that are of 
higher order in the momenta we can express x m terms of the source: 



X 



2ttT' 



(2.19) 
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Similarly we can verify that the following expression is correct up to higher order 
terms in the momenta: 



diag(2\/2i + iq - 4(1 - y/2£)u, 2y/2i -4q- 4(1 - y/2t)uj)x (2.20) 



Inverting the matrix and expanding on small momenta we can express x i n terms 
of the source: 

X = diag(-t^- - 4g + 4(1 - V2t)uj, -i^- + Aq + 4(1 - V2l)u)tp (2.21) 

2.3 Computation of the diffusion constant 

Given that, for the Rarita-Schwinger field, the on-shell action vanishes , the only 
contribution to the correlator will come from the boundary action: 

Sbdy = M j d 4 xV^hh ij ^ l ¥ (2.22) 

where h is the induced metric on the boundary and *I f = i^T 1 . Performing the 
field redefinitions <j> = 7 2 V>2 + 7 3, 03 an d tj = ^ 2 4>2 — 7 3, 03 an d setting cj) = 0, ipi = 0, 
ip4 = since we are only interested in the transverse part we get: 



Fourier transforming and taking the limit u = e — > we have: 

r] ] _{k)ri + {k) + rjUk)ri-{k) (2.24) 



Shdv = — — r— / d 4 k\J —h 
y 2 ttTR J 

The correlator that we want to calculate is not hermitian and therefore we should 
keep only the first term in the action to extract the correlator. This procedure is 
equivalent to the recipe proposed in [12]. Substituting the expansion of the field on 
the boundary the action becomes: 

S bdy = y (vrTfl) 3 J d*k (f^ + e Vy x ) ( 2 .25) 

To substract the divergencies in the limit e — > we introduce the counterterm: 

Set = St% = ~y^ TR ) 3 j d*te -1/ Vx (2-26) 
The counterterm is written in covariant form as: 

S c t = -^ I d 4 xV^hfj-Ari- (2-27) 
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UJ)<p 



where f/_ = rfT 1 and $y l = T^d^. This counterterm does not introduce any finite 
correction to the correlator since there is no u 7 / 4 term in the boundary expansion 
of ?7_. The correlator is calculated from the renormalized action. 

Sren — S^dy S& — 

~ (ttTR) 3 J diag(-^ - 4q + 4(1 - + 4g + 4(1 - yftl) 

The retarded correlator of the operator dual to r\ is therefore: 

G n = y(vrr J R) 3 diag(-^-4 (? + 4(l-\/2^)a;,-^ +4^ + 4(1 -y/2l)u) (2.28) 

where we have absorbed the extra e power in a field redefinition. The normalization 
of the action is calculated in the appendix and is found to be: 

N 2 

N = -f (2.29) 

From the analysis of the projections of the correlators we can deduce that with i 
indices summed • = 4G?: where GfH- is the correlator that comes in the Kubo 



formula (A. 11) for the diffusion constant. The energy density is: 

3tt 2 



. N 2 T 4 

8 c 

We finally compute the diffusion constant from the Kubo formula and find: 

D s = (2.30) 
This result is in agreement with the diffusion constant calculated in [7]. 

2.4 Running of the correlator 

In this section we give an alternative derivation in the spirit of the membrane 
paradigm approach to the question of universality of the diffusion constant [6] . The 
idea is that the correlator can be formally computed on any radial slice of the 
geometry. When evaluated on the boundary it corresponds to the real field-theory 
correlator, and one can derive a flow equation that expresses the running from the 
boundary to the horizon. In favorable cases, at low frequency the flow is trivial and 
the UV correlator is equal to the value at the horizon. 



Let us consider again the Rarita-Schwinger equation (2.14); at zero momentum 
and frequency, it becomes 

ri' = {-A + B) v+ , 7]'_ = {-A - S)r?_ (2.31) 
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where B 



and m = 3/2 for the gravitino. The ratio w = r/ + /r]^ satisfies 



0, 



w 



2Bw, and we have the flow equation 



T ' = du ( e -2fB(u)duV+\ =Q 



V- 



(2.32) 



The integral converges at the horizon, so we have J-(u = 1) = r]+/rj-(u = 1) = —i. 
At the boundary we have J- ~ u^ m rj + /r]^ = x/^i Du t the coefficient of proportion- 
ality requires the computation of J," ^/g U u so we cannot express it purely in terms 
of horizon data. We find 



D, 



9vrT 



exp 



-m log e — 2m 




(2.33) 



This formula reproduces the result (2.30), and allows to compute the diffusion con- 



stant in any asymptotically AdS black hole background. However, as we will see, in 
presence of a finite charge density for a U(l) under which the gravitino is charged 
(typically an R-symmetry), the equations of motion are more complicated and do 
not become diagonal even at zero momentum. 



3 Supercharge diffusion at finite density 
3.1 The background 

We want to consider the change in the diffusion constant when there is a finite 
background charge density. In N=4 SYM there is an SU(4) R-symmetry group 
under which the supercurrents transform in the fundamental representation. It is 
possible to give a vev simultaneously to three charges in the Cartan subalgebra of 
SU(4). The corresponding solution is an asymptotically AdS black hole carrying 
electric charges, known as the STU black hole. It was first constructed as a solution 
of d = 5,M = 2 gauged supergravity with 2 matter multiplets. The construction of 
the theory is explained in [13] and the solutions of the bosonic part are presented 
in [T3] . The solution can be embedded in type IIB supergravity, the uplift to 10 
dimensions is given in [15] but for our purposes it is sufficient to work with the 
5-dimensional solution. 

The theory contains three abelian gaueg fields Aj and three scalar fields subject to 
the constraint X 1 X 2 X S = 1 , with the metric on the scalar manifold given by: 

G IJ = 2 diag[(X 1 ) 2 , (X 2 ) 2 , (A 3 ) 2 ] (3.1) 

and the scalar potential is V = 2^2j 1/X 1 . The metric of the black hole with 
translationally invariant horizon is: 

ds 2 = _ n -2/3(^f 2 nl/3 (KTLf 2 2 ^ n l/3_tf du 2 
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f(u) = U{u)-u 2 \\{l + n I ) H i (u) = 1 + kiu 



Kl 



1=1 



where qj are the three charges. In the zero charge limit this metric is exactly (2.8). 
The scalar and gauge fields are given by: 

xI = irr% A * = * TLu irr\ II ^ + «j) ( 3 - 2 ) 

Hj(u) -H/(n) 
The Hawking temperature of the background is: 

Th _ 2 + KX + K2 + K3- K 1 K 2 K 3 ^ 
2^/(1 + Kl)(l + K2)(l + K 3 ) 

The energy density, entropy density, charge density and chemical potentials are 
respectively: 

3ir 2 N 2 T 4 



Kj 



TT 2 N 2 T 3 J J , ., 

s = — 2 — na+Kj) 1 / 2 



1=1 

3 



(3.4) 



1=1 

pi=— §— v^n( i+ ^) 1 
j=i 

1 + KI j=i 

Imposing the condition of thermodynamic stability on the background implies a 
restriction on the charges of the background [16J 

2 — K\ — K2 — ^3 + K1K2K3 > (3.5) 

The relevant part of the Lagrangian for the gravitino is [13J: 

1 

-/ 
2 

(3.6) 

where I = 1,2,3, i, a = 1,2. hf, V a and V are functions of the scalar fields , 4> l are 
the scalars fields on the constrained scalar manifold and /f is the vielbein on this 
manifold. The covariant derivative acts on a spinor as 



where Vi are constants and should be equal to 1/3 for the Hj to be properly nor- 
malized [H]. V is given by: 

Vo = \g U yj (3.7) 

and F^ u = d^A^ — d u A^ are the gauge field strengths. The equations of motion for 
the gravitino are: 



+ ^T^uVo - ^X a V a = 



The supersymmetry transformations act on the gravitino as: 



(3.8) 



5V, = 2V + TO, - ^r")^ p e + ^ E V ( 3 ' 9 ) 



where 



= -I — * i 

and we neglected the terms containing the spinors A in the two previous formulas. 
We are again interested in the transverse component of the gravitino, which can 
be defined as r\ = T x ^ x — T y ^ y . Using the same strategy as in the zero-charge 
case, one can show that the transverse part of the gravitino is decoupled, both 
from the longitudinal part and from the spinors A. The equation for the transverse 
component is therefore: 



with 



+ F{u) - ^V&"WVr4 + -^i^JHFit + ^v&'Vo) V = 

(3.10) 



tP- ■ % 1 ,-4 

W = —i><u -ji v + iqi 
f(u) = l9 k + 3 sL 

4 gtt 4 g xx 



(3.11) 



Note that the transverse component of the gravitino is gauge invariant and that we 
can again split our equation in two different systems by projecting on the eigenspaces 

Of 723- 



3.2 Analytic solution for small charge 

The system of differential equations obtained above cannot be solved analytically. 
Since calculating the diffusion constant requires knowing the retarded correlator for 
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zero momenta we set the momentum to zero and attempt to find a solution as an 
expansion in the charges. We consider in turn the cases where only one charge is 
turned on, and when there are three equal charges. 



One charge 

In this case we take K\ = K, k,i = K3 = and attempt to find a solution of the form: 
rj(u, k) = r/ (0) (u) + (u) + k r/ 2) (u) + . . . (3.12) 

Solving the system of equations near the horizon we identify the exponents: 



tf ± = 4±l^±M±^ (3.13) 
± 4 4 2+k y ' 

The choice of U- yields incoming- wave boundary conditions. Near the boundary 
we find the following behavior: 

m = ^(u, n)u^ + X (u, kW" + *( w , k) V^+^) + ^ 3 u 7/4 lQg u + 

77 3 = <f>(u, K)u 1/4: + 4>(i0, K)u 5/i + ... 

This expansion implies that there will be a logarithmic divergence of the correlator 
that would need to be regularized with an appropriate counterterm. However this 
will not affect the diffusion constant that is computed from the imaginary part of 
the correlator, so we will ignore this divergence. 

The integration constants that appear at every order in the yfH expansion (3.12) 
can be fixed at zeroth order by imposing incoming boundary conditions at the 
horizon and at the i-th order by imposing 77® = on the horizon. The near 
boundary expansion up to order ^fn for the two components i^j 

771 = n 7 / 4 [2" 3 / 4 i + v^(a + (3 log u)\ 
m = n 1 / 4 (2 3 / 4 - iy^ 7 ) 

where a, /3, 7 are real constants. Notice that there is no term going like u 3//4 so we 
don't need to regularize the boundary action. In terms of the source <f> and the 
normalizable mode x the correlator is then simply given by: 

= —^-4> = — i+ ^-4^ + 2Wi l0gU J^ + - (3 - 16) 

The coefficient of y/H is real and therefore the diffusion constant is not corrected to 
the order ^[k. 



2 For q=0 the system for (771,773) is identical to the system for (772,774) 
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Three charges 



We consider here the case of three equal charges, namely k± = K2 = K3 = K. For 



this special choice of the charges the thermodynamic stability condition (3.5) is 



saturated for k = 2 but it is not violated even for k > 2. The Hawking temperature 



(3.3) on the other hand becomes negative for k > 2 and so in order to have a 
physical theory we must restrict the value of the charge to k < 2. For k = 2 we 
have an extremal black hole. 

We attempt to find a solution as an expansion in ^J~k as in the case of one charge. 
The indices near the horizon are 

4 4(2-k)VTT^ 

and the incoming wave boundary condition corresponds to V- . The near boundary 
expansion is written in this case: 



// \ 3/4 , / \ 7/4 , // \ 3v k(1 + k) — 8w 7 / 4 
?7l = U)(f){UJ, K)U 1 +XV UJ ,K)U' +(p(Ul,K) — - u 1 Iogn+... 

8 

7? 3 = 4>(U, H)U 1/A + 4>(UJ, K)n 5/4 + . . . 

(3.18) 

The calculation done for the one charged case goes through for the three charges 



case and the end result is similar to (3.15) with different numerical values for the 



constants a, (3, 7. In this case too, then, the diffusion constant does not receive 
corrections to order ^[k. 

The explicit computation of higher orders in this expansion becomes too compli- 
cated, therefore we revert to a numerical analysis. 



3.3 Numerical Solution 

In order to solve numerically the system of equations we need to provide initial 
conditions in form of a near horizon expansion where we have imposed the incoming 
boundary conditions. We will work with the components (771, 773) since the equations 
for the other two are exactly identical for q = 0. We solve the system near the 
horizon and find a expansion of the form: 

Vhor(u, It, K) = 7]^ r (LJ, K) + y/l-U ^(w, K) + ... (3.19) 

Using the expansion above as an initial condition we integrate the system of equa- 
tions from the horizon to the boundary. Then we can extract the correlator from 
the near boundary region of the solution. 

The result for the one-charge case is presented in figure [T] At k = it reproduces 
correctly the value calculated analytically for zero charge. 
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Figure 1: The diffusion constant as a function of the charge for the one charge 
case K\ = K, Kl = K3 = . 

In figure [2] we present the diffusion constant as a function of T//i and find a 
behavior similar to the one observed in |17| for an AdS-RN solution of N=2, D=4 
supergravity. For large values of T j \x the diffusion constant approaches the zero 
density value which is represented by the horizontal line. 
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Figure 2: The diffusion constant as a function ofT/fi for the one charge case. 
The horizontal line represents the zero density analytical result. 



In figure [3] we give the result for three charges as a function of k, and in figure [4] 
as function of T//j>. We observe that the behavior is qualitatively similar to that of 
the one charge case for high temperature, but the diffusion constant (or rather the 
combination TD S ) approaches a finite limit as the temperature goes to zero. 
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Figure 3: The diffusion constant in the three equal charges case as a function of 
the charge. 
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Figure 4 : The diffusion constant as a function ofT/jj, for the three charges case. 
The horizontal line represents the zero density analytical result. 



4 Conclusions 

We have considered the hydrodynamic limit of the correlators of supercharges in 
N=4 SYM at zero and finite R-charge density. From the transverse part of the cor- 
relator the supercharge diffusion constant can be extracted using a Kubo formula. 
In the zero charge case we could confirm the result found previously from the cal- 
culation of the supersymmetric sound attenuation, using the longitudinal channel 
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of the correlator. 

The main motivation for the present work has been to answer the question of the 
presence of universality in the hydrodynamic supersymmetric sector of theories with 
a holographic dual. The conclusion appears to be negative: the diffusion constant 
depends on the value of the U(l) charge density. We have not been able to establish 
the precise form of the dependence, in particular there does not seem to be a simple 
relation between the diffusion constant and the thermodynamic quantities. However 
qualitatively we found a behavior very similar to the one found in the 4-dimensional 
AdS-RN black hole. 

The absence of universality is perhaps not surprising, since if we consider the form 
of the equation (3.10) we see that the transverse gravitino is coupled to the scalars 
and gauge fields that are running in the solution. 

We observed that the dimensionless diffusion constant TD S tends to a finite value 
for the extremal black hole. It would be interesting to know if this is true in all 
cases of extremal black hole, and perhaps find if this value can be predicted by some 
sort of attractor mechanism. 
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A Kubo formula derivation 

The constitutive relation [9] which relates the density of charge to the spatial com- 
ponents of the supercurrent 

S l a = -D B &p a + D a (^d jP ) a - ~( 7 Vp) a (A.l) 

Combining with the continuity equation dtp a + diS* a and Fourier transforming in 
space yields: 

[(d t + D s k 2 )6? - ic s {h?irf)£] P fi(t, k) = (A.2) 

where c s = P/e. The solution of this equation in the limit of small k is: 

p a (t, k) = e~ D ° kH [bi cos(kc s t) + i(k j 7jiy a sm(kc s t)]pp(0, k) (A.3) 

From this solution we can deduce the following supercharge density correlators in 
the small k limit since C a Q,(0,0) is fixed by the supersymmetry algebra. 

Ca&(t,k) = f d s xe- i% ' s {p a (x,t)pa(d,0)) = 

J „. (A ' 4) 

= -£e~ DskH [(i'y 1 ) adi cos(kc s t) + (k j jj) a a sin(kc s t)} 

3 When we write 7^ we mean the upper right block of the 7 matrix. We made the replacement 
a — > ij 1 and o % — > —i'f with respect to [?J]. 
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We define the following correlator: 

C a a{t,x) = -i9{t)( Pa {t,x)pa(0,d)) 

We can now calculate the Fourier transform: 
Ca&(w,k) = -i J dt J d z xe^ t - it H{t){ Pa (t,x)p 6l {{)^)) 
= -i I dte iut 6{t)C a a{t,k) 



(A.5) 



iD s k 2 + u 



(c s k) 2 + (D s k 2 - ioj) 



r(«7 1 )ad + 



(c s k) 2 + (D s k 2 - ioj) 2 



The retarded correlator of the supercurrents is defined as: 

CL(t,*) = -i8(t)({p a (t,x),p A (0,d)}) 
From the spectral decomposition of C and C R it follows: 

Im C a a(uj, k) = ~Yl e~ fiEn (n\p a (0, k)\m) {m\p&(0, 0)\n)5{uj + E n - E m ) 



(A.6) 
(A.7) 



lmC^(oj,k) 



7T 

z 



Therefore we get: 



-l3E n 



ImC^(co,k) = (l + e- puJ )ImC a a(oj,k) 



[l + e~^){n\p a (0,k)\m){m\pa(0,0)\n)5(u + E n -E„ 



(Ai 



Taking the small k and then the small oj limit of (A.6) and using equation (A. 8) we 
get: 

1 



eD s = - lim 
4 w-+o 



lim^(i, ')"" hn('Hu.k) 



(A.9) 



Using the continuity equation we can deduce a relation between the supercharge 
and the supercurrent correlators, namely: 



u 2 C^(oj,k) = k 2 kkG l L{u,k) = 



(A.10) 



where G l ^(t,x) = — i0(t)({Sa(t,x), S 3 & (0, 0)}) and in the last equality the i indices 
are summed. Using this relation we can write the Kubo formula in terms of the 
supercurrent retarded correlator: 



1 



lim 



12 uj^o 



lim^ 1 )"" JmG"Juj,k) 



(A.ll) 
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B Normalization of the boundary action 



In this appendix we fix the normalization of the boundary action using the zero 
temperature calculation of the super stress energy tensor correlator. Exceptionally 
in this appendix we use the conventions of |18j . 

The leading term of the super energy-momentum tensor OPE has the following 
structure up to a constant A [19]: 



"(s 2 s 2 ) 

where 

i 



S a 

Ja a (z)J j3/ j(0) = A . „ „. 2 (B.l) 



•S 



aa 



(x - x') aa + -[o a (§ - e% + e' a (e- e'u 



s aa = {x- x') aa + % -[e a (e - 6% + e' a (e- e') a ] 

The super energy-momentum tensor is a superfield containing the axial current 
(at zeroth order in 9) , the supercurrent (at linear order in 9) and the energy- 
momentum tensor (at second order in 9). To fix the proportionality constant we 
can calculate the energy-momentum tensor correlator by applying twice the operator 
D a ] + [Dp, Dp] J a a) on the super energy-momentum tensor OPE and 
compare with the energy momentum correlator |20j: 

T^uTpv = — j- p CT ( 3fiK&v\ + ^iX^uk — -SuvSkx) (B.2) 

7T \X \ \ _■ / 



where J" = 6" — 2x ili x v j\x\ 2 . From the comparison we deduce that: 

A = -- i Nj (B.3) 



7T J 



Knowing the exact expression fot the super energy-momentum tensor OPE we can 
compute the supercurrent correlator by applying the operator D^D'^ and doing the 
proper symmetrizations. Then we can compare with the zero temperature correlator 
of the supercurrent calculated in [10] : 



16AT TTpC j° <i x A / x p x v 



S«S« = ~^K^T [v P u ~ 2^ j (B.4) 

where Il£ = (Sfc — j7/i7 p ) and find the value of the normalization of the action: 

N 2 

A/-=% (B.5) 
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